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ABSTRACT 

Some massive stars possess strong magnetic fields that confine plasma in the circumstellar environment. 
These magnetospheres have been studied spectroscopically, photometrically and, more recently, interferomet- 
rically. Here we report on the first firm detection of a magnetosphere in continuum linear polarization, as a 
result of monitoring of a Ori E at the Pico dos Dias Observatory. The non-zero intrinsic polarization indicates 
an asymmetric structure, whose minor elongation axis is oriented 150.°0 east of the celestial north. A modula- 
tion of the polarization was observed, with a period of half of the rotation period, which supports the theoretical 
prediction of the presence of two diametrally opposed, co-rotating blobs of gas. A phase lag of -0.085 was de- 
tected between the polarization minimum and the primary minimum of the light curve, suggestive of a complex 
shape of the plasma clouds. We present a preliminary analysis of the data with the Rigidly Rotating Magneto- 
sphere model, which could not reproduce simultaneously the photometric and polarimetric data. A toy model 
comprising two spherical co-rotating blobs joined by a thin disk proved more successful in reproducing the 
polarization modulation. With this model we were able to determine that the total scattering mass of the thin 
disk is similar to the mass of the blobs {2M\,/Mi = 1 .2) and that the blobs are rotating counterclockwise on the 
plane of the sky. This result shows that polarimetry can provide a diagnostic of the geometry of clouds, which 
will serve as an important constraint for improving the Rigidly Rotating Magnetosphere model. 
Subject headings: circumstellar matter — stars: individual (Sigma Ori E) — stars: magnetic field 



1. INTRODUCTION 

The helium-strong star crOrionisE (HD 37479; B2Vpe; 
my = 6.66) has long been known to possess a circum- 
stellar magnetosphere that is formed by wind plasma that 
is trapped by a strong dipolar magnetic field (wlO kG, 
iLandstreet & Borra Il978l) . The presence of this magneto- 
sphere can be inferred by the eclipse-like dimmings on its 
light curve, which is thought to occur when p lasma clouds go 
in front of the star twice every rotation cycle (iTownsend et al.l 
llOOl iHesseretaniTOTTt fGroote & Hungedll982h . This star 
has been the testbed for important advancements in our the- 
oretical understanding of these magnetospheres. The Rigidly 
Rotating Magnetosphere (RRM) models of this star were suc- 
cessful in reproducing the o bserved variability in emission 
line profiles and photometry (iTownsend et al.ll2005l) . 

Polarization is a very useful technique that allows one to 
probe the geometry of the circum stellar scattering rnateri al 
without angularly resolving it (e.g.. lBrown & McLean 1977). 
iKemp & Herman! ([1977) carried out the first polarimetric ob- 
servations of the a Ori system, but their results where largely 
inconclusive because of the small values of the intrinsic polar- 
ization. In this paper we present the results of high-precision 
(cr ^ 0.01%) polarization monitoring of crOriE, that has re- 
sulted in the first firm detection of the polarization modulation 
produced by a co-rotating magnetosphere. 

2. OBSERVATIONS 
Electronic address: |carciofi@usp.br| 



Broad-band linear polarization data were obtained from 
August, 2010 to September, 2011, using the lAGPOL im- 
age polarimeter attached to the 0.6-m Boiler & Chivens 
telescope at Pico dos Dias Observatory (OPD/LNA). We 
us ed a CCD camer a with a polarimetric module described 
in iMagalhaes et alj (Il996l) . consisting of a rotating half- 
waveplate and a calcite prism placed in the telescope beam. 
In each observing run at least one polarized standard star was 
observed in order to calibrate the observed position angle. 
HD 23512 and HD 187929 were used as polarized standards. 

Although typical polarimetric observations consists of 8 
consecutive wave plate positi ons (hereafter WPP) separated 
by 22.°5 (C arciofi et al.ll2007l) . in this work we chose to use 
16 consecutive WPPs in order to increase the signal-to-noise 
ratio. Also, to reach the high accuracy needed to monitor the 
modulation of the polarization, a large number of frames (typ- 
ically from 20 to 100) were obtained at each WPP. Since the 
exposure times ranged from 0.3 to 1 s per frame, the tempo- 
ral resolution of our observations vary from 7 to 34 min, in- 
cluding telescope and instrument overheads. These values are 
short enough to temporally resolve the polarimetric variation 
over the 1 . 19 d period of the star. Detail s of the data reduction 
can be found in lMagalhaes et alj (11984 . 

3. INTERSTELLAR POLARIZATION 

iKemp & HermanI d 19771) report similar values of the polar- 
ization for (J Ori AB, a Ori C and a Ori D, which suggests that 
all these stars belong to the same physical system, and that the 
measured polarization is of interstellar (IS) origin. 
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Fig. 1. — Variation of the intrinsic polarization of a Ori E compared with photometric variations in the y filter. Top panels: y filter photometric data from 
IHesser et alj {1977). folded in phase. Middle panels: intrinsic P (red), Q (black) and U (blue). Bottom panels: intrinsic position angle. Left plot: The polarization 
data is displayed in the equatorial reference frame. The polarization of a Ori AB was used as a measure of the IS polarization. Middle plot: Same as the left 
plot, rotated such that = 0. Ri^ht plot: Polarization compared to the single-scattering model. The solid line is for our best-fitting model with = 70° and 

the dashed line is for i= 110°. For this plot we used the estimate for the IS polarization and d made in Sect. [6] 



We measured the BVRI polarization of cr Ori AB and 
C and our results are largely consistent with the ones 
of iKemp & Hermaiil (IT977). We used our measured V- 
band polarization of aOri AB (P = 0.351(15)%; 86.7°; Q = 
-.348(15)%; U = 0.040(15)%) as an estimate of the IS polar- 
ization towards a Ori E. 

4. RESULTS 

The V-band cr Ori E observations are shown in the left 
plot of Fig. [T] The data was folded i n phase (0) using the 
ephemerides of iTownsend et alJ ( l2010l) . Shown is the intrin- 
sic polarization, calculated by subtracting the IS polarization, 
estimated using a Ori AB, from the observed Q and U Stokes 
parameters. The polarization has a double-peaked structure, 
with maxima of about 0.07% occurring at phases around 0.3 
and 0.8 and minima occurring around phases 0.1 and 0.6. We 
note that th e values reported here are in broad agreement with 
the data of iKemp & HermanI jl977l) . but a more quantitative 
comparison is hampered by the insufficient accuracy of the 
latter data. 

There is a clear anti-correlation between the photomet- 
ric and polarimetric variations, with the polarization minima 
roughly coinciding with the photometric maxima and vice- 
versa, but there are important differences between the two 
curves. One such difference lies in their symmetry: the pho- 
tometry is clearly asymmetric, with the secondary minimum 
occurring at about phase 0.4 (the primary minimum occurs 
at = by definition), and the behavior of the two inter- 
minima being very different from each other The polarization 
curve, however, seems to be roughly symmetric with respect 
to w 0.6. That the asymmetry is seen in the photometry, but 
not in the polariza tion, is further supporti ng evidence that - 
as hypothesized in iTownsend et alj (l2005h - the asymmetric 



light curve arises due to contamination of the circumstellar 
signal by a photospheric signal. That is, the photometric vari- 
ations come from a combination of magnetospheric eclipses 
and photospheric abundance inhomogeneities, and therefore 
lack rotational symmetry; but conversely, because the polar- 
ization variations come solely from magnetospheric electron 
scattering, we do see rotational symmetry in their case. 

In order to determine the location of the polarization 
minimum, we fitted the data between (j> = 0.45 and 0.75 
with a para^bola, using the Levenberg-Marquardt algorithm 
dMarquarddl 1 963h . We find that the minimum of the intrinsic 
polarization occurs at ^min = 0.61 ±0.03. To determine the 
location of the maximum, the same procedure was applied to 
the data between = 0.7 and 1. According to this fit the max- 
imum of the polarization occurs at (p^^j^ = 0.84 ± 0.04. 

To further assess the symmetry of the data we fitted the po- 
larization data with the following function 

P((j)) = Po+Acos[An{(j)-S)], (1) 

which implicitly assumes a period of half of the orbital pe- 
riod for the polarimetric data. The result of the fit gives Pq = 
0.0471 ±0.0009%,A = 0.021 ±0.001 % and 5 = -0.17±0.01, 
with a reduced = 0.67. The fitted function is displayed as 
the black line in the middle plot of Fig. [T] 

The data seems to be well represented by Eq. ([1), which 
indicates that, contrary to the photometry, the polarization 
curve is roughly symmetric. From Eq. dO, the position of 
the minima are (j) = 0.08 and 0.58, and the maxima are at 
(j) = 0.33 and 0.83. Those values are in good agreement with 
the phases derived above via a parabola fitting. The mini- 
mum and maximum values of the polarization are, respec- 
tively, Po-A = 0.026 ±0.001% andPo+A = 0.068 ± 0.001%. 

An important feature of the intrinsic polarization of a Ori E 
is that it is never zero. This indicates that there is some degree 
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Fig. 2. — Modeling of the intrinsic polarization of a Ori E using the RRM 
model (observations are in red). The only free parameter is the maximum 
number density in the magnetosphere, which was set to lO'^cm"^ (sohd 
hnes) to reproduce the depth of the eclipses and 2.5 X lO" cm"' (dotted hnes) 
to reproduce the amplitude of the linear polarization. 

of asymmetry of the scattering material (as seen projected 
onto the plane of the sky) throughout the entire rotation pe- 
riod. 

The data shown in the left plot of Fig. [T] are in the equatorial 
frame. The weighted-average of all the data gives a value of 
(2int) = 0.023% and (i/im) = -0.040%, which, in turn, results 
in an average position angle of (6'int) = 150.°0. This value tells 
us the (average) direction of the minor axis of the asymmet- 
ric structure. Furthermore, the small variation of the position 
angle (Fig.[Tl left) tells us that the direction of the minimum 
elongation axis changes little (« 10-15°, at most) as the star 
rotates. 

This last point is better understood if we rotate the intrinsic 
polarization by (6'int) so that the average value of U is zero 
(Fig. [T] right plot). In this new frame, we see that the am- 
plitude of the rotational modulation of U is small. If we as- 
sume that the symmetry axis is parallel to the rotation axis, 
we conclude that the projected axis of symmetry of the mag- 
netosphere varies little as the star rotates. As discussed below, 
this has quite important consequences for the RRM model. 

Another noteworthy feature of the polarization curve is that 
there is a phase shift of -0.085 between the primary mini- 
mum of the lig htcurve and the polarization minimum. The 
RRM model of Townsend et alj j2005h predicts the existence 
of two plasma clouds co-rotating with the star If the clouds 
were symmetrical and the photospheric flux were homoge- 
neous, there should be no phase difference between the po- 
larization and photometric curves; the observed phase shift, 
therefore, is likely a result of cloud asymmetries and photo - 
spheric abundance inhomogeneities. 

5. RRM MODEL 

iTownsend et alj ( l2005h apphed the RRM model to cr Ori E 
with good success. This model predicts the accumulation of 
circumstellar plasma in two co-rotating clouds, situated in 
magnetohydrostatic equilibrium at the intersections between 
the magnetic and rotational equators. Comparison with the 
available data (photometry and Ha line profiles) showed that 
the model provided a good quantitative description of the cir- 
cumstellar environment. 

In an attempt to reproduce the observed polarization mod- 
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Fig. 3. — Geometric conception of the "dumbbell + disk" model, to scale. 

ulation of cr Ori E, we fed the predicted density d istribu- 
tion of the RRM model of Townsend et al. (2005) to the 
HDUST radiative transfer code (Carciofi & Bjorkman 2006). 
Since all stellar and g eometrical parameters were taken from 
iTownsend et alJ ( l2005i) . the only free parameter in the model 
is the density scale of the magnetosphere. 

The results of this modeling are shown in Fig.|2] We could 
not find a model that matched simultaneously both the pho- 
tometry and the polarimetry. A higher density model (solid 
line) that matches the depth of the photometric eclipses pre- 
dicts a polarization amplitude that is three times larger than 
what is observed. Conversely, a lower density model that re- 
produces the amplitude of the polarization fails to reproduce 
the photometric amplitude. 

Furthermore, the detailed shape of the polarization curve, 
in particular the position angle variation, is not well matched. 
Since the polarization postion angle is sensitive to the geom- 
etry of the scattering material, this discrepancy indicates that 
the primary difficulty with the basic RMM model is the shape 
of its density distribution. 

6. SINGLE SCATTERING MODEL 

We now explore what changes must be made to the basic 
RRM model to better reproduce the shape of the polarimet- 
ric variability. The RRM predicts the existence of diametrally 
opposed plasma clouds confined in the rotational equator and 
a diffuse plasma disk close to the magnetic equator Based on 
this idea, we developed a simple toy model for a Ori E that 
consists of two elements: 1) a thin uniform disk tilted by an 
angle if) from the rotational equator, and 2) a pair of spher- 
ical blobs situated in the equatorial plane at the intersection 
between this plane and the magnetic equator, thus forming a 
"dumbbell-like" structure. This configuration is illustrated in 
Fig [3] for several orbital phases and for a viewing angle of 
70°. 

Since the polarization is small, we may perform the cal- 
culation of the Q and U Stokes parameters as a function of 
orbital phase for a given inclination angle / following a sim- 
ple single-scattering approach (valid only in the optically thin 
limit). The scattered flux is first deter mined in the frame of the 
star using the formalism described in iBiorkman & BiorkmanI 
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( 1994 Eg. 6) and then rotate d to the frame of the observer 
(iBjorkman & Bjorkmanlll994 Eq. 20). 

For this initial study, we adopted 5 fixed parameters, listed 
in Table [T| wh ose values were based on the models of 
iTownsend et al.l ('2005). They are the stellar radius, R^,, the 
distance of the blobs from the star, db, the size of the blobs, R^, 
the geometrical thickness of the disk, H^, and the inclination 
angle, /. The adopted value for /?biob is somewhat arbitrary, 
since in the limit of /?biob ^ R* and in the single-scattering 
approximation, a large, tenuous blob is roughly equivalent to 
a dense, small blob. For the same reason, the disk geometrical 
thickness is also somewhat arbitrary: if the single-scattering 
approximation holds, a thinner, denser disk is equivalent to a 
thicker, less dense disk. In other words, what really controls 
the polarization level of each component is their total scat- 
tering mass. For the inclination angle we adopt two possible 
values: / = 70°, corresponding to a projected counterclock- 
wise rotation of the blobs on the sky and 180° -70° = 110°, 
corresponding to a clockwise rotation. 

In this analysis we made no previous assumptions about 
the IS polarization; i.e., we fitted directly the observed Q and 
U Stokes parameters. The adopted modeling procedure is as 
follows: 

1 . Choose values for electron number density in the blob 
and disk and n^), the phase lag, 5, between the pho- 
tometric primary minimum (0 = 0) and the orbital phase 
where the blobs are aligned with the line of sight, the 
tilt angle of the disk, ijj, and the position angle of the 
equatorial plane in the plane of the sky, 0. These 5 pa- 
rameters define a unique model for the intrinsic Q and 
U in the equatorial frame. 

2. Choose values for Qis and Uys and use them to calculate 
synthetic observed Q and U parameters. 

3. Calculate the reduced of the model. 

4. Repeat steps 1 and 3 for several hundred of values for 
the 7 free parameters and choose the model with mini- 
mum x^. 

The parameters of our best-fitting model (reduced = 2.8) 
are listed in Table [T] The estimated errors are for a 95% con- 
fidence interval. A notable result is that the IS polarization, 
independently obtained by this analysis, is consistent (within 
the errors) with the estimate made using a Ori AB as a mea- 
sure of the IS. Both the orientation of the equatorial plane 
in the sky {9) and the phase lag {5) are consistent with the 
estimate made in Sect. |4] The resulting blob density corre- 
sponds to an electron optical depth of r = 0.13, which does 
not violate the optically thin assumption. It is worth men- 
tioning that our estimate for the electron density in the blobs, 
lO'^e" cm"-^, is in good a gree ment with the values derive d by 
iGroote & Hunged (1 1982b and ISmith & Bohlended (l2007l) us- 
ing the Inglis-Teller formula. From the best-fitting param- 
eters one can calculate the mass of each model component. 
Assuming a molecular weight of 0.6, which corresponds to a 
fully ionized gas with solar chemical composition, the mass 
of each blob is Mb = 6.0 x 1O"'^M0 and the mass of the disk 
is Md = 1.0 X 10"" M0. Therefore, we obtain that both com- 
ponents have similar masses (2Mb/Md = 1.2). 

In Fig.[Tl right plot, we show the observations corrected by 
the IS polarization of Table [T] The solid curves correspond 
to the best-fitting single scattering model. The two minima 



TABLE 1 
Parameters of the single 
scattering model 



r cQ'3rn c tcr 


V 1 

aue 


lype 


/ (degree) 


70 or 110 


fixed 


R* (Rq) 


4.3 


fixed 


4 (R*) 


2.4 


fixed 


Rb (R*) 


1/3 


fixed 


Hi (R.) 


0.01 


fixed 


{e cm"') 


1.0^-6 X 10'2 


fitted 


(e~ cm"') 


2.7±1 X 10'2 


fitted 


'0 (degree) 


28 


fitted 


eis(%) 


-0.35 ±0.01 


fitted 


Uis (%) 


0.025 ±0.010 


fitted 


9 (degree) 


150 ±7 


fitted 


<5 


-0.17 ±0.02 


fitted 



of the model polarization occur at the phases for which the 
two blobs are alined with the line of sight. At these phases, 
the model Q parameter is minimum and the U parameter is 
maximum (negative for « (5 and positive for (p « 0.5 + (5), 
because the disk dominates the polarized flux at these phases. 
Conversely, the polarized flux is dominated by the blobs at 
the phases at which the blobs are on the side of the star {(j> « 
0.25 + S and positive for 4> ~ 0.75 + 6). 

An interesting result is that the data allowed us to firmly 
establish that the blobs are rotating counterclockwise on the 
plane of the sky. In the right plot of Fig. [T] the solid curves 
show the best-fitting model for / = 70° (reduced = 2.8) 
while the dotted curves correspond to the best fitting model 
with/= 110° (reduced x^ = 5.1). The former is clearly a much 
better fit to the data. 

An alternate comparison between the data and the single- 
scattering model is made in Fig.|4] where we plot the observed 
(i.e., uncorrected for IS polarization) U vs. 2 parameters. The 
blue points correspond to the phase range S-(S + 0.25), the 
orange ones from ((5 + 0.25)-((5 + 0.50), the green ones from 
((5 + 0.50)-((5 + 0.75) and the red ones from (S + 0.75)-S. In 
the QU plane the behavior of the data is quite different from 
the first half-cycle to the second: while in the first half (blue 
and orange points) the data forms a nearly-circular outer loop, 
in the second half (green and red points) there is an inward 
incursion. The model reproduces this behavior 

It is important to point out that the two-component model 
(blob-ndisk) adopted by us is required to produce such a be- 
havior in the QU plane. For instance, a model with only the 
blobs would form a symmetric double ellipse in the QU plane, 
while a model with only the tilted disk would result in a shape 
resembling a crescent moon. The observed curve, composed 
of an outer loop and an inner incursion, can only be obtained 
by a combination of these two components. 

7. CONCLUSIONS 

The polarimetric data presented here for a Ori E has an ac- 
curacy of at least one order magnitude better than those shown 
in the work of Kemp & Herman (1977) and thus represent 
the first firm detection of the polarization modulation of a co- 
rotating magnetosphere. The intrinsic polarization, calculated 
using a Ori AB as a measure of the IS polarization, displays 
a roughly symmetric sinusoidal curve ranging from 0.02 to 
0.08%, with a period of half of the orbital period and a phase 
lag of -0.085 between the polarization minimum and the pri- 
mary minimum of the light curve. The polarization position 
angle is also variable. Its mean value of 150.°0 indicates that 
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However, the data was well reproduced by a two- 
component model, consisting of two spherical, diametrally 
opposed blobs and a disk tilted with respect to the rota- 
tional equator. This model was physically motivated by 
the RRM model, which predicts accumulation of centrifu- 
gally supported wind plasma in a disk-like structure with the 
highest densities (corresponding to the blobs) at the points 
where the magnetic equator intersects the rotational equator 
This model provided useful constraints on the total scatter- 
ing mass of each component. Also, it provided an indepen- 
dent estimate for the IS polarization, as well as for the ori- 
entation of the equatorial direction in the plane of the sky. 
Our best-fitting model suggests that the mass of each blob 
is Mb = 0.6 X 1O"'^M0 and the mass of the disk is Md = 
1.0 X 1O~"M0, which indicates that both components have 
similar masses (2Mb/Md = 1 .2). Also, the tilt angle of the disk 
with respect to the plane of the sky is 28 deg, simil ar to the an- 
gle be tween the rotational and magnetic axes ( Townsend et alj 
l2005h . Finally, the model allowed us to determine that the 
blobs are rotating counterclockwise on the plane of the sky. 



Fig. 4. — QU plot of the observed data (points with error bars). The orbital 
phase correspondence of the points is as follows: blue, (5-((5+0.25), orange, 
(5 + 0.25H<5 + 0.50). green, (5 + 0.50H(5 + 0.75), and red, (<5 + 0.75)-<5. The 
filled circles show the best-fitting model of Table [T] The IS polarization of 
crOriAB is indicated. The grey diamond shows the IS of Table[T] 

the minor elongation axis of the asymmetric structure around 
a Ori E is oriented 150.°0 east of the celestial north. 

A preliminary analysis made with the RRM model 
jrbwnsend et al.ll2005h suggests that it is difficult to fit both 
the photometry and the polarization simultaneously. This 
seems to imply that the geometry of the plasma clouds, as 
predicted by the current version of the model, is incorrect. 
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